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Abstract 

The effect of the presence of a finite impedance surface on the wave propagation properties of 
a two-dimensional periodic array of rigid cylinders with their axes perpendicular to the surface is 
both numerically and experimentally analyzed in this work. In this realistic situation both the 
incident and the scattered waves interact with these two elements, the surface and the array. The 
interaction between the excess attenuation effect, due to the destructive interference produced by 
the superposition of the incident wave and the reflected one by the surface, and the bandgap, due 
to the periodicity of the array, is fundamental for the design of devices to control the transmission 
of waves based on periodic arrays. The most obvious application is perhaps the design of Sonic 
Crystals Noise Barriers. Two different finite impedance surfaces have been analyzed in the work 
in order to observe the dependence of the wave propagation properties on the impedance of the 
surface. 

PACS numbers: 43.20.+g, 43.35. +d 
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I. INTRODUCTION 



Research on the transmission properties of Phononic Crystals (PC),-'^ defined as periodic 
lattices of elastic scatterers embedded in an elastic medium with different physical properties, 
have attracted increasing interest in last decade. Exhaustive studies have been performed 
in order to know the underlying physics in their main control properties of sonic and elastic 
waves. Refraction," invisibility,- bandgaps,-"^ or recently, control of phase properties,- 
are some topics related with the exploitation of the wave propagation properties of PC. 
In these years, several theoretical methods have been developed to analyze the physical 
properties of these systems. Plane Wave Expansion method (PWE)^ allow to find the bands 
structure of the PC transforming the wave equation into an eigenvalue problem that can be 
solved for each Bloch vector k, in the irreducible first Brillouin zone. Extended Plane Wave 
Expansion (EPWE)^ii^ is a general method for analyzing the complex dispersion relation 
of PC. Multiple Scattering Theory (MST)i^"— is a self consistent method to estimate the 
acoustic pressure in a point taking into account the incident pressure on the PC and the 
scattering due to all the scatterers. Finally, Finite Elements Method (FEM)^ can be used 
to characterize periodic media where the geometries of the scatterers or the existence of 
some new effects different than scattering are introduced in the problem. 

In a real situation, scatterers should be hold on a finite impedance surface and the 
refiections of both, the incident and the scattered waves on this surface have to be taken 
into account when the transmision properties of PC are analyzed. However, this effect has 
not been considered in the previous models. As a first approximation to this problem, recent 
worksi^'^ present both semi analytical and numerical models to analyze the infiuence of a 
finite impedance surface on the transmission properties of arrays of rigid cylindrical scatterers 
embedded in air, usually called sonic crystals (SC), with their axes parallels to the surface 
(horizontal scatterers). This situation has the most tractable geometry to develop a semi 
analytic model. However, as we have mentioned, the real case involves arrays of cylinders 
with their axes perpendicular to the surface which means a three dimensional (3D) geometry 
problem. This case implies to assume the complete geometry of the composite resulting in 
very complex analytical developments. On the other hand, numerical simulations of these 
3D geometries would imply a long computational time. 

The goal of this work is to present numerical predictions and experimental results to 
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FIG. 1. (Color online) Experimental set up used 



show the variation of the transmission properties of SC devices when the cyhnders are per- 
pendicular to the finite impedance surface. To do this we develop a simplified 3D numerical 
model based on FEM and we check the results with accurate experimental results obtained 
in controlled conditions. The interaction between the excess attenuation effect, due to the 
destructive interference produced by the superposition of the incident wave and the wave 
reflected by the surface, and the bandgap, due to the periodicity of the array, has been found 
of fundamental interest for the design of devices based on periodic arrays devoted to control 
the transmission of waves. The study is focused on the airborne transmission case. The 
results obtained here allow to enhance the knowledge about the sonic behaviour of these 
systems in a more realistic way. 

The work is organized as follows: First of all we explain the used experimental set up and 
the developed numerical model. Then a discussion about the obtained numerical results and 
a comparison with the experimental ones is shown. Finally, the concluding remarks show 
the main conclusions of the work. 

II. EXPERIMENTAL SET UP AND NUMERICAL ANALYSIS 
A. Experimental set up 

The experimental results have been obtained under controlled conditions in an echo-free 
chamber. In fig. [1] one can observe the experimental set up used. We have constructed a 
SC made of 20 (5x4) acoustically rigid cylinders with radius r = 0.09 m and 1.20 m length, 
arranged in a square array with lattice constant a = 0.22 m. With this geometry the first 
pseudogap at FX direction appears between 500 and 980 Hz, as we will explain later. 



Two types of surfaces have been considered. The first one is a 0.03 m thick wooden board, 
theoretically considered as an acoustic rigid surface. The second one is a 0.1 m foam sheet, 
considered as a porous material with finite impedance characterized using a two-parameter 
impedance modeU^ with flow resistivity cr=250 kPa s/m^ and porosity at the surface a = 1 
m~^. To experimentally study the interaction between the SC and the surface, we have 
taken measurements on a vertical line of points with length y = 1 m (starting at 0.05 m 
from the surface) behind to the SC, as one can see in fig. [H We have measured this line 
in steps of Ay = 0.02 m in order to analyze in-depth the interaction between the surfaces 
and the bandgaps of SC by varying the height of the microphone. The distance between the 
vertical line of measurements and the center of the nearest row of cylinders is c = 0.5 m, 
and the horizontal distance between the sound source and the measurement line is b = 2.51 
m. A source emitting white noise is placed at d = 0.3 m distance from the ground. 

In our experimental set up, the acquisition is done by a prepolarized free-field microphone 
1/2" Type 4189 B&K. The microphone position is varied by using a Cartesian robot which 
controls the movement along the three axes (OX, OY and OZ), installed in the ceiling of 
the echo-free chamber. The robot is called by us 3DReAMS (3D Robotized e-Acoustic 
Measurement System), and it has been designed to sweep the microphone through a 3D grid 
of measuring points located at any trajectory inside the chamber. When both the acoustic 
source and the robotized system are turned off, the microphone acquires the temporal signal. 
This signal is saved on the computer and then, using the Fast Fourier Transform (FFT), 
one can obtain the power spectra, the frequency response or the sound level measurement. 
National Instruments cards PCl-4474 and PCI- 7334 have been used together with the Sound 
and Vibration Toolkit and the Order Analysis Toolkit for Lab VIEW for both the data 
acquisition and the motion of the robot. 

B. Numerical analysis: Finite Element Method 

To numerically solve the problem we consider the geometry shown in fig. A line 

source is centred at point ro = {x,y). The emitted acoustic field by this source at point r 
is po = Ho{k\f — fo\), where Hq{x) is the Hankel function of the first kind and zero order. 
A row of 4 scatterers separated by the lattice constant of the array, a, is enclosed between 
two completely reflecting walls which are parallel to the axis of the scatterers. These walls 
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FIG. 2. (Color online) (a) Numerical solution domain, (b) Numerical simulation for the scattering 
problem with completely rigid ground at 1000 Hz. The real value of the total pressure is plotted, 
Re{P). 

are also separated by a lattice constant, a. It is worth noting that the incident wave is not 
reflected by these walls, however the scattered waves from the cylinders are reflected by the 
walls reproducing the effect of a semi infinite periodic slab of 4 rows of scatterers. Then 
a reduced volume is used here to reproduce the effect of a semi infinite structure. This 
procedure is only valid for plane waves and line sources emitting cylindrical waves. 

The commercial software COMSOL 3. 5. a is used in this work to obtain the numerical 
predictions. The discrete domain where the solution is calculated is surrounded by Perfectly 
Matched Layers (PML)i^ which are an efficient alternative for emulating the Sommerfeld 
radiation condition in the numerical solution the scattering problems. The finite impedance 
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surface where the cyhnders are hold is modelled using the impedance boundary condition. 
We solve the problem for the scattered waves as the incident wave ps is itself solution of 
the wave equation. Thus, the total acoustic field is obtained adding the incident wave to 
the scattered one, then the PML only acts on the scattered wave. The solved problem has 
10^ degrees of freedom and 83.7 10^ elements have been used to solve it. An example of the 
numerical solution at 1000 Hz is shown n fig. |2]d. 

The situation without the SC is very well known: the total field can be characterized by 
both the incident and the reflected waves. This last one can be obtained from the reflection 
coefficient of the finite impedance surface and the image method. Then, without the SC, 
the total acoustic field is^^ 



Ptot=PO+Ps (1) 

where po is the incident wave and ps = R{ro, r; iy)Ho{k\r' — r'o|) represents the reflected one 
being R{fQ, r; u) the reflection coefficient of the finite impedance surface. On the other hand, 
the parameter used in this work to characterize the attenuation properties of the complete 
system (SC and finite impedance surface) is the Insertion Loss (IL) defined here as the 
difference between the sound level recorded with and without the SC considering always the 
presence of the finite impedance surface. Then, 

/L = 201og^ (2) 

where P is the total acoustic field with both the surface and the SC. We would like to remark 
that we have taken as reference for the IL the situation where the surface is present. 

Sometimes the dispersion relation of the periodic system is useful to compare the effect 
of the finite size on the scattering problems. Thus, the bands structures of the designed SC 
have been also calculated using the PWE. 961 plane waves have been used in the calculation 
as this number produced a good convergence of the method for the case analysed here. We 
fix our attention in the range of frequencies included in the first pseudogap due to the fact 
that this region has important relevance in the field of the acoustic barriers which is perhaps 
the most exploited application of SC for audible noise in airborne transmission. 
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III. RESULTS AND DISCUSSION 



First of all we analyze the case of the finite impedance surface without the SC. Fig. [3] 
shows both numerical predictions (1) and experimental results (2) of the dependence of the 
excess attenuation (EA) on the height of the measurement point and the frequency for the 
rigid (a) and absorbent (b) surfaces. In the EA regions a destructive interference between 
the direct and the reflected waves on the surfaces is produced, and as a consequence an 
attenuation region is observed. Note that the position of these regions depends on both the 
height and on the frequency. 
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FIG. 3. (Color online) EA Theoretical analysis in free field conditions for a rigid surface (al) 
and finite impedance surface (bl). EA Experimental measurements performed with 3DReAMS 
for the two previous cases: rigid surface (a2) and finite impedance surface (b2). The geometrical 
conditions of both the theoretical and the experimental results have been defined in Section [TI] 

On the other hand, in fig. [3] one can check the fairly good agreement between the numeri- 
cal and the experimental results, as a proof of the accurate precision of our experimental set 
up. The existing differences in the maps for the case of rigid surface (wood) can be explained 
in terms of the assumption of acoustically rigid surface in the numerical predictions. In the 
case of the finite impedance surface, EA appears at lower values of both frequencies and 
heights than in the case of rigid ground. 

Once the EA effect of the surface is characterized we analyze the interaction between 
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FIG. 4. (Color online) IL theoretical analysis of the SC-rigid surface interaction (al) and SC- 
finite impedance surface interaction (bl). Experimental results obtained with SDReAMS for the 
two previous cases: SC-rigid surface (a2) and SC-finite impedance surface (b2). The results are 
obtained with the geometrical conditions explained in Section II. B. The bands structure has been 
calculated using the PWE method, with the next parameters: Pcj/«mders=2700 kg/m^, /3aj,.=1.3 
kg/m^; Ccyiinders =6400 m/s and Cair=340 m/s, and with 961 plane waves. We have focused our 
attention in the properties of the first band gap at FX, between 500 and 980 Hz 



SC and both the rigid and the absorbent surfaces. Figs. H^l and 11^2 show the numerical 
predictions and the experimental results respectively for the case of the acoustically rigid 
surface. In this figures one can see that the maximum values of the IL appear inside the 
range of frequencies corresponding to the first band gap at FX direction, represented in the 
figures with the vertical dotted lines as a prolongation of the edges of the bandgap shown in 
the band structures calculated using PWE method. For this case the interaction produces 
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a positive effect on the EA of the surface, that means, an increasing of the attenuation 
properties of the surface. On the other hand, there are some regions on the map where 
the IL is negative for both several heights and frequencies, which mean a reinforcement of 
the pressure level compared to the case of rigid surface alone which is considered as the 
reference. Note that this reinforcement is produced by the combined effect of the SC and 
the surface, and it means that the attenuation level is lower than in the case of the rigid 
surface alone, although some attenuation still exists compared with the free field conditions. 

Similar effects can be seen in figs. Hbl and |Dd2 for the case of finite impedance surface. 
The range of frequencies where the IL is higher is included in the limits of the band gap 
at rX direction, marked with vertical dashed lines. Moreover, some regions in the map 
show a reinforcement of the acoustical level, although here the size of these areas is lower 
than in the rigid surface case. All these results are in good agreement with the theoretical 
results showed for the SC-surfaces interaction in the two-dimensional case (horizontal 

cylinders) . 

Up to now, we have compared the dependence of the EA on the different types of surfaces 
when a SC is introduced. In figure 5 we show the opposite analysis that explains the influence 
of the surfaces in the attenuation properties of SC. In this figure, one can compare the IL 
spectra taking as reference the free field or the corresponding surface for three heights of 
the microphone and for the two surfaces analyzed. The dashed lines in fig. [5] show the IL 
calculated considering ptot in eq. [2] the acoustic field considering only the source and the 
finite impedance surface. The continuous lines represent the IL calculated considering ptot 
the corresponding for the free field (without the surface). In the case of rigid surface (figs. 
[5^1, [5]d1 and[5tl) one can see the increasing of the attenuation properties of the SC when 
the rigid surface is considered. In the case of the finite impedance surface there is not a 
relevant variation in the size of the FX pseudogap when the surface is considered. 

IV. CONCLUDING REMARKS 

We have shown both numerically and experimentally the interaction between the attenu- 
ation properties of both a SC and a finite impedance surface. We have focused our analysis 
on the case where the axis of the rigid cylinders that form the SC are perpendicular to the 
surfaces. For this realistic case we have found a good agreement between both the numerical 
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FIG. 5. Experimental IL spectra for three heights of the microphone (y=91 cm (a); y=65 cm (b) 
and y=31 cm (c)). Continuous and dashed hues represent the IL of SC without groud and with 
ground respectively. Figure marked with 1 (2) represents the measurements for rigid (absorbent) 
ground. One can see in the case of the rigid surface its effect on the attenuation properties of SC 
(marked with arrows). In the case of a finite impedance surface the existence of the surface does 
not affect the SC attenuation properties. 



predictions and the experimental data measured using SDReAMS. Both the EA, due to the 
presence of the finite impedance surface, and the band gap of the SC interact in order to 
show a combined effect that depends on the frequency and on the relative position between 
the measuring point and the position of the source. We have analyzed this interaction under 
the point of view of the attenuation properties of the SC, finding that for the case of rigid 
surfaces the combined effect is relevant and it has to be taken into account. However, for 
finite impedance surfaces with high absorption values, the interaction practically disappears 
and the attenuation properties of the SC in the measuring points are only dependent on the 
periodicity of the array as well as on the number of scatterers. The results shown in this 
work can be useful to design attenuation devices based on periodic materials as SC acoustic 
barriers as well as for the design of experiments in which the effect of the surface could 
interact with the transmission properties of SC. 
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